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We have performed a thorough analysis of the constraints which can be put on neutrino parameters 
from cosmological observations, most notably those from the WMAP satellite and the 2dF galaxy 
survey. For this data we find an upper limit on the sum of active neutrino mass eigenstates of 
^2 m v < 1.0 eV (95% conf.), but this limit is dependent on priors. We find that the WMAP and 
2dF data alone cannot rule out the evidence from neutrinoless double beta decay reported by the 
Heidelberg- Moscow experiment. In terms of the relativistic energy density in neutrinos or other 
weakly interacting species we find, in units of the equivalent number of neutrino species, iV„, that 
N„ = 4.0±!;? (95 % conf.). When BBN constraints are added, the bound on N v is 2.6±%\j (95 % 
conf.), suggesting that N v could possibly be lower than the standard model value of 3. This can 
for instance be the case in models with very low reheating temperature and incomplete neutrino 
thermalization. Conversely, if N v is fixed to 3 then the data from WMAP and 2dFGRS predicts that 
0.2458 < Yp < 0.2471 (95% conf.), which is significantly higher than the observationally measured 
value. The limit on relativistic energy density changes when a small v e chemical potential is present 
during BBN. In this case the upper bound on N v from WMAP, 2dFGRS and BBN is N„ < 6.5. 
Finally, we find that a non-zero ^ m„ can be compensated by an increase in N v . One result of this 
is that the LSND result is not yet ruled out by cosmological observations. 

PACS numbers: 14.60.Pq,95.35.+d,98.80.-k 



I. INTRODUCTION 

Neutrinos exist in equilibrium with the electromag- 
netic plasma in the early universe, until a temperature 
of a few MeV. At this point the weak interactions freeze 
out and neutrinos decouple from the plasma. Shortly af- 
ter this event, the temperature of the plasma falls below 
the electron mass, and electrons and positrons annihi- 
late, dumping their entropy into the photon gas. This 
heats the photon gas while having no effect on neutri- 
nos, and the end result is that the photon temperature 
is larger than the neutrino temperature by the factor 
(11/4) 1 / 3 ~ 1.40. Since the present day photon tem- 
perature has been measured with great accuracy to be 
2.728 K, the neutrino temperature is known to be 1.95 
K, or about 2 x 1(T 4 eV. Since the heaviest neutrino has 
a mass of at least about 0.04 eV it must at present be 
extremely non-relativistic and therefore acts as dark mat- 
ter. The contribution of a single neutrino species of mass 
rn v to the present day matter density can be written as 
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so that for a neutrino mass of about 30 eV, neutrinos 
will make up all of the dark matter. However, this 
would have disastrous consequences for structure for- 
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mation in the universe, because neutrinos of eV mass 
have very large free streaming lengths and would erase 
structure in the neutrino density on scales smaller than 
if s ~ 1 m~l Y Gpc completely. This leads to an over- 
all suppression of matter fluctuations at small scales, an 
effect which is potentially observable. 



A. Absolute value of neutrino masses 

The absolute value of neutrino masses are very dif- 
ficult to measure experimentally. On the other hand, 
mass differences between neutrino mass eigenstates, 
(mi, TO2, 77^3 ) , can be measured in neutrino oscillation ex- 
periments. Observations of atmospheric neutrinos sug- 
gest a squared mass difference of Sm 2 ~ 3 x 10~ 3 eV 2 
While there are still several viable solutions to 
the solar neutrino problem from solar neutrino observa- 
tions alone, the large mixing angle (LMA) solution gives 
by far the best fit with Sm 2 ~ 5 x 10~ 5 eV 2 HEt Re- 
cently the KamLAND reactor neutrino experiment has 
announced a positive detection of neutrino oscillations 
indicating that the LMA solution is indeed correct [24| . 

In the simplest case where neutrino masses are hi- 
erarchical these results suggest that mi ~ 0, 777,2 ~ 
t>TO so iar, and rri3 ~ 6m a tmos P heric- If the hierarchy is 
inverted H H H IH 0| one instead finds m 3 ~ 0, 

771 2 ~ ^atmospheric, and 777-1 ~ ^atmospheric- HoW- 

ever, it is al so p ossible that neutrino masses are de- 
generate [13, Ea E3, EH EE E3, EM E3, HH, 

mi ~ m-2 ~ m-3 3> <5m a t m ospheric, hi which case oscillation 
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experiments are not useful for determining the absolute 
mass scale. 

Experiments which rely on kinematical effects of the 
neutrino mass offer the strongest probe of this overall 
mass scale. Tritium decay measurements have been able 
to put an upper limit on the electron neutrino mass 
of 2.2 eV (95% conf.) |2^|. However, cosmology at 
present yields an even stronger limit which is also based 
on the kinematics of neutrino mass. As discussed be- 
fore any structure in the neutrino density below the free- 
streaming scale is erased and therefore the presence of 
a non-zero neutrino mass suppresses the matter power 
spectrum at small scales relative to large scale, roughly 

by ap/p ~ -8n v /n m M- 

This power spectrum suppression allows for a determi- 
nation of the neutrino mass from measurements of the 
matter power spectrum on large scales, as well as the 
spectrum of CMB fluctuations. This matter spectrum is 
related to the galaxy correlation spectrum measured in 
large scale structure (LSS) surveys via the bias parame- 
ter, b 2 (k) = P g (k)/P m (k). Such analyses have been per- 
formed several times before |78l l79j , most recently using 
data from the 2dFGRS galaxy survey IHiSlp. These 
investigations found mass limits of 1.5-3 eV, depending 
on assumptions about the cosmological parameter space. 

In a seminal paper it was calculated by Eisenstcin, 
Hu and Tegmark that future CMB and LSS experiments 
could push the bound on the sum of neutrino masses 
down to about 0.3 eV [7j|. The prospects for an abso- 
lute neutrino mass determination was discussed in fur- 
ther detail in Ref. |2!j where it was found that in fact 
the upper bound could be pushed to 0.12 eV (95% conf.) 
using data from the Sloan Digital Sky Survey and the 
upcoming Planck satellite. 

More recently the new WMAP data, in conjunction 
with large scale structure data from 2dFGRS has been 
used to put an upper bound on the sum of all neutrino 
species of £ m v < 0.7 eV (95% conf.) H. 

However, the exact value of this upper bound de- 
pends strongly on priors on other cosmological param- 
eters, mainly Hq. In the present paper we calculate the 
upper bound on m v from present cosmological data, 
with an emphasis on studying how the bound depends 
on the data set chosen. 

In addition to their contribution to the cosmological 
mass density neutrinos also contribute to the cosmolog- 
ical energy density around the epoch of recombination. 
Neutrinos which have mass smaller than roughly 3T rec , 
where T rec ~ 0.3 eV is the temperature of recombina- 
tion, will act as fully relativistic particles when it comes 
to CMB and large scale structure. 

In the standard model there are three light neutrino 
species with this property. However, these particles are 
not necessarily in an equilibrium Fermi-Dirac distribu- 
tion with zero chemical potential. It is well known that 
the universe contains a non-zero baryon asymmetry of 
the order rj — nB n — ~ 10~ 10 . A neutrino asymmetry of 
similar magnitude would have no impact on cosmology 



during CMB and LSS formation, but since the neutrino 
asymmetry is not directly observable it could in princi- 
ple be much larger than the baryon asymmetry. Such a 
neutrino asymmetry would effectively show up as extra 
relativistic energy density in the CMB and LSS power 
spectra. 

Another possibility for extra relativistic energy is that 
there are additional light species beyond the standard 
model which have decoupled early (such as the graviton 
or the gravitino). 

From the perspective of late time evolution at T <■ 1 
MeV it is customary to parametrize any such additional 
energy density in terms of N v [3(j, the equivalent num- 
ber of neutrino species. In Section III we discuss bounds 
on N v from the present WMAP and 2dFGRS data, com- 
bined with additional information on other cosmological 
parameters from the Hubble HST key project and the 
Supernova Cosmology Project. 

However, as will be discussed later, a non-zero neu- 
trino chemical potential can have an effect on big bang 
nucleosynthesis which is profoundly different from simple 
relativistic energy density if it is located in the electron 
neutrino sector. This possibility will be further discussed 
in Section III. 

Another important point is that any entropy produc- 
tion which takes place after BBN, but prior to CMB for- 
mation will only be detectable via CMB and LSS ob- 
servations. One such example is the decay of a hypo- 
thetical long-lived massive particle at temperatures be- 
low roughly 0.01 MeV. 

The general outline of the paper is as follows: In Sec- 
tion II we discuss the likelihood analysis as well as the 
data sets used in the analysis. In Section III we present 
the numerical results of the analysis in terms of an up- 
per bound on Y m v, a confidence interval on N u , and 
a bound on the neutrino chemical potentials. Finally, 
section IV contains a discussion and conclusion. 



II. LIKELIHOOD ANALYSIS 

The extraction of cosmological parameters from cos- 
mological data is a difficult process since for both CMB 
and LSS the power spectra depend on a plethora of differ- 
ent parameters. Furthermore, since the CMB and matter 
power spectra depend on many different parameters one 
might worry that an analysis which is too restricted in 
parameter space could give spuriously strong limits on a 
given parameter. We discuss this point in further detail 
below. 

For calculating the theoretical CMB and matter power 
spectra we use the publicly available CMBFAST package 
|31j. As the set of cosmological parameters we choose 
Q, m , the matter density, the curvature parameter, f^, the 
baryon density, Hq, the Hubble parameter, n s , the scalar 
spectral index of the primordial fluctuation spectrum, r, 
the optical depth to reionization, Q, the normalization 
of the CMB power spectrum, b, the bias parameter, and 
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finally the two parameters related to neutrino physics, 
Q, y h 2 and N u . We restrict the analysis to geometrically 
flat models, i.e. f2 = £l m + £!a = 1. 

In principle one might include even more parameters 
in the analysis, such as r, the tensor to scalar ratio 
of primordial fluctuations. However, r is most likely 
so close to zero that only future high precision exper- 
iments may be able to measure it. The same is true 
for other additional parameters. Small deviations from 
slow-roll during inflation can show up as a logarithmic 
correction to the simple power-law spectrum predicted 
by slow-roll 13^13^1341 or additional relativistic energy 
density [HlISWM El El E| could be 
present. However, there is no evidence of any such effect 
in the present data and therefore we restrict the analy- 
sis to the "minimal" standard cosmological model plus 
neutrino physics. 

In this full numerical likelihood analysis we use the 
free parameters discussed above with certain priors deter- 
mined from cosmological observations other than CMB 
and LSS. In flat models the matter density is restricted 
by observations of Type la supernovae to be Q m = 
0.28 ± 0.14 HH, and the HST Hubble key project has 
obtained a limit on H of 72 ± 8 kms" 1 Mpc -1 pr| . 

We calculate bounds on neutrino physics, both with 
and without, the constraints on f2 m and Hq. In the cases 
where we do not use the SNI-a and HST priors we use 
simple top-hat priors instead, 0.1 < O m < 0.5 and 0.5 < 
h < 0.85. 

Particularly the prior on Hq is of great importance in 
constraining neutrino parameter space. 

The numerical marginalization over parameters other 
than Q, v h 2 was performed using a simulated annealing 
procedure ■ A full description of how x 2 is calculated 
can for instance be found in lial. 



A. LSS data 

At present, by far the largest survey available is the 
2dFGRS [H of which about 147,000 galaxies have so 
far been analyzed. Tegmark, Hamilton and Xu have 
calculated a power spectrum, P(fc), from this data, which 
we use in the present work. The 2dFGRS data extends 
to very small scales where there are large effects of non- 
linearity. Since we only calculate linear power spectra, we 
use (in accordance with standard procedure) only data 
on scales larger than k = 0.2h M pc" 1 , where effects of 
non-linearity should be minimal 29]. Making this cut 
reduces the number of power spectrum data points to 
18. 



B. fc-dependent bias 

In all present parameter estimation analyses it is as- 
sumed that the bias parameter, b 2 (k) = P g (k)/ P m (k), is 
independent of the scale, k. 



However, many independent simulations find that bias 
is in fact quite strongly scale dependent 0, |S2 in the 
non-linear regime. In the linear regime bias is expected to 
be constant, and the asymptotic value = Hmfc_>o&(fc) 
is reached as a scale of roughly k ~ 0.1 — 0.2h Mpc -1 . 
This means that at scales larger than fc cut bias should be 
very close to scale-independent, and that we can therefore 
use a single parameter, b, to describe it. 

C. CMB data 

The CMB temperature fluctuations are conveniently 
described in terms of the spherical harmonics power spec- 
trum 

Ci = (\a lm \ 2 ), (2) 

where 

AT 

—(e,<i>) = J2aimYi m (e,<t>). (3) 

Im 

Since Thomson scattering polarizes light there are ad- 
ditional powerspectra coming from the polarization 
anisotropics. The polarization can be divided into a curl- 
free (E) and a curl (B) component, yielding four indepen- 
dent power spectra: Cr,z, Ce.i, Cb,i and the temperature 
i?-polarization cross-correlation Cte,i- 

The WMAP experiment have reported data only on 
C T ,i and C TE ,u as described in Ref. [H |H |H |H |H3 

We have performed the likelihood analysis using the 
prescription given by the WMAP collaboration which in- 
cludes the correlation between different C/ 's [2^, IBIj. I54I 
I55ll56j |. Foreground contamination has already been sub- 
tracted from their published data. 

In parts of the data analysis we also add other CMB 
data from the compilation by Wang et al. |74[ which 
includes data at high I. Altogether this data set has 28 
data points. 

III. NUMERICAL RESULTS 
A. Neutrino masses 

We have calculated x 2 as a function of neutrino mass 
while marginalizing over all other cosmological parame- 
ters. This has been done using the data sets described 
above. In the first case we have calculated the constraint 
using the WMAP C T ,i combined with the 2dFGRS data, 
and in the second case we have added the polarization 
measurement from WMAP. Finally we have added the 
additional constraint from the HST key project and the 
Supernova Cosmology Project. It should also be noted 
that when constraining the neutrino mass it has in all 
cases been assumed that N v is equal to the standard 
model value of 3.04. In section III.D we relax this condi- 
tion in order to study the LSND bound. 
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The result is shown in Fig. 1. As can be seen from 
the figure the 95% confidence upper limit on the sum of 
neutrino masses is m » < 1-01 e V (95% conf.) using the 
case with priors. This value is completely consistent with 
what is found in Ref. |82j where simple Gaussian priors 
from WMAP were added to the 2dFGRS data analysis. 
For the three cases studied we find the following limits: 



Y,m v < 1.01 eV 
J2m v < 1.20 eV 
Y,m v < 2.12 eV 



However, it is somewhat higher than the upper limit 
of Y, m v < 0.7 eV found in the WMAP analysis |28] . 
There are several reasons for this: First, we do not use 
Ly-a forest data in our analysis. While adding them to 
the analysis clearly would improve the limit the interpre- 
tation of Ly-a data is still somewhat controversial |84j . 
Therefore the limit derived here can be viewed as more 
robust, albeit also more conservative. The second rea- 
son is that we use a completely free bias parameter. It 
well known that a determination of the neutrino mass 
is to some extent degenerate with the bias parameter b 
so that a precise determination of b would shrink the al- 
lowed range even further. However, the bias parameter is 
probably the most poorly understood parameter in the 
analysis of large scale structure [74| and therefore im- 
posing a particular prior on b could yield too restrictive 
results. This is the reason why we do not show any result 
for \ 2 with a 6-prior. 

Also, for accurate CMB and LSS data sets, the main 
degeneracy is not with the bias parameter, but rather 
with the Hubble parameter. This was also noted in 
Ref. [i^. 

In the middle panel of Fig. 1 we show the best fit value 
of Hq for a given fl^h 2 . It is clear that an increasing value 
of m u can be compensated by a decrease in Hq. Even 
though the data yields a strong constraint on Sl m h 2 there 
is no independent constraint on f2 m in itself. Therefore, 
an decreasing Hq leads to an increasing £l m . This can be 
seen in the bottom panel of Fig. 1. 

When the HST prior on Ho is relaxed a higher value 
of ^2 m v is allowed, in the case with only WMAP and 
2dFGRS data the upper bound is Sl v h 2 < 0.023 (95% 
conf.), corresponding to a neutrino mass of 0.71 eV for 
each of the three neutrinos. 

This effect was also found by Elgar0y and Lahav |82j in 
their analysis of the effects of priors on the determination 

Of Y^ m v 

However, as can also be seen from the figure, the addi- 
tion of high-? CMB data from the Want et al. compila- 
tion also shrinks the allowed range of ^2 m v significantly. 
The reason is that there is a significant overlap of the 
scales probed by high-Z CMB experiments and the 2dF- 
GRS survey. Therefore, even though we use bias as a free 
fitting parameter, it is strongly constrained by the fact 
that the CMB and 2dFGRS data essentially cover much 
of the same range in fc-space. 
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FIG. 1: The top panel shows \ 2 as a function of 5^m„ for 
different choices of priors. The dotted line is for WMAP + 
2dFGRS data alone, the dashed line is with the additional 
Wang et al. data. The full line is for additional HST and SNI- 
a priors as discussed in the text. The horizontal lines show 
A\ 2 = 1 and 4 respectively. The middle panel shows the best 
fit values of Ho for a given ^ m„ . The horizontal lines show 
the HST key project la limit of H = 72 ± 8 kmsMpc" 1 . 
Finally, the lower panel shows best fit values of Q m - In this 
case the horizontal line corresponds to the SNI-a la upper 
limit of Q m < 0.42. 



It should be noted that Elgar0y and Lahav [g^] find 
that bias does not play any role in determining the bound 
on m v ■ At first this seems to contradict the discussion 
here, and also what was found from a Fisher matrix anal- 
ysis in Ref. ji^. The reason is that in Ref. [g^, redshift 
distortions are included in the 2dFGRS data analysis. 
Given a constraint on the amplitude of fluctuations from 
CMB data, and a constraint on £l m h 2 , this effectively 
constrains the bias parameter. Therefore adding a fur- 
ther constraint on bias in their analysis does not change 
the results. 

Neutrinoless double beta decay 

Recently it was claimed that the Heidelberg-Moscow 
experiment yields positive evidence for neutrinoless dou- 
ble beta decay. Such experiments probe the 'effective 
electron neutrino mass m ee — | U 2 ^m Vj \. Given the 



5 



uncertainties in the involved nuclear matrix elements 
the Heidelberg-Moscow result leads to a mass of m ee — 
0.3 — 1.4 eV. If this is true then the mass eigenstates are 
necessarily degenerate, and ~ 3m ee . Taking the 

WMAP result of m v < 0.70 eV at face value seems to 
be inconsistent with the Heidelberg-Moscow result |8fij . 
However, already if Ly-a forest data and a constraint on 
the bias parameter is not used in the analysis the upper 
bound of J2m„ < 1.01 eV is still consistent. For this 
reason it is probably premature to rule out the claimed 
evidence for neutrinoless double beta decay. 



B. Neutrino relativistic energy density 

For the case of the effective number of neutrino 
species we have in all cases calculated constraints in the 
{Q,f } h 2 ,N l/ ) plane, while marginalizing over all other pa- 
rameters. The reason for this is that for Big Bang Nu- 
cleosynthesis purposes these are exactly the important 
parameters. Therefore, to combine CMB, LSS, and BBN 
constraints the marginalization over flbh 2 should not be 
performed. Furthermore, when constraining N v we have 
always assumed that m v ~ so that the neutrino mass 
has no influence on cosmology. 

We start out by investigating constraints on N v from 
CMB and LSS data alone. In Fig. 2 we show A% 2 
for a global fit to N v which marginalizes over all other 
cosmological parameters. The overall best fit for the 
WMAP T and TE data, combined with the Wang et 
al. compilation, the 2dFGRS data, the HST key project 
data on Hq, as well as the SNI-a data on f2 m , has 
X 2 = 1467.6 for a total of 1395 degrees of freedom. This 
gives x 2 /d.o.f. = 1.05 which is entirely compatible with 
the best fit WMAP value for the standard ACDM model 
of x 2 /d.o.f. = 1.066. We also show constraints for two 
other analyses. The first is for WMAP and 2dFGRS data 
alone and the second for WMAP data alone. The bounds 
for the three cases are 
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FIG. 2: x 2 as a function of ^ N v for different choices of priors. 
The dotted line is for WMAP data alone, the dashed line is 
with the additional Wang et al. and 2dFGRS data. The full 
line is for additional HST and SNI-a priors as discussed in the 
text. The horizontal lines show A% 2 = 1 and 4 respectively. 
The middle panel shows the best fit values of Ho for a given 
N v . The horizontal lines show the HST key project lcr limit of 
Ho — 72 ± 8 kmsMpc -1 . Finally, the lower panel shows best 
fit values of f2 m . In this case the horizontal line corresponds 
to the SNI-a la upper limit of Q m < 0.42. 
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These bound are entirely compatible with those found by 
Crotty, Lesgourgues and Pastor |73|. and much tighter 
than the pre- WMAP constraints. 

The constraints derived here are also compatible with 
what is found by Pierpaoli [T^ , where are assumption of 
spatial flatness was relaxed. 

In the lower panels of Fig. 2 we show the best fit values 
of Hq and f2 m for a given value of N u . The main point to 
note is that the constraint on N v is strongly dependent 
on Hq. This was also found in Ref. [H- With only CMB 
data and the weak top-hat prior on Hq the bound on N v 
is very weak. Adding the HST Key Project prior on Hq 
cut away a significant amount of parameter space both 
at low and high N v . Adding the 2dFGRS and Wang et 
al. data mainly has the effect of shifting the best fit value 



to higher N v , but also cuts away the low N v values, an 
effect also seen in Ref. |7^. 

In Fig. 3 we show constraints on (f2b/i 2 , Af„) for the 
full data set described above. The best fit value for f^/i 2 
is 0.0233, which is equivalent to the value found in the 
WMAP data analysis. In the 2-dimensional plots the 
68% and 95% regions are formally defined by A\ 2 = 2.30 
and 6.17 respectively. Note that this means that the 68% 
and 95% contours are not necessarily equivalent to the 
same confidence level for single parameter estimates. 

It should be noted here that in addition to an upper 
bound on N u there is also a 3. Oct confidence detection 
of N v > 0. This is in concordance with the pre- WMAP 
data from which a non-trivial lower bound on N v could 
also be derived. 

Adding BBN information - In the case where all the 
relativistic energy density contained in N v is produced 
prior to BBN, a BBN constraint can be added to the 
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FIG. 3: 68% and 95% confidence contours in the (Qbh 2 , N„) 
plane for the WMAP TT and TE data, combined with the 
2dFGRS data, the HST data on H and the SNI-a data on 



FIG. 4: 68% and 95% confidence contours in the {Q. b h 2 ,N v ) 
plane for the same data sets as in fig. 3, but with the addition 
of BBN data. The lined contours are the 68% and 95% regions 
for BBN data alone. 



CMB and LSS constraint without any problems. In prac- 
tice we have used abundances of He-4 and D to make 
constraints in the (Sibft 2 , iV„) plane. We use the follow- 
ing values for the primordial abundances |63l l64j | 

D/H = 2.78±g;|| x 1(T 5 (4) 
Y P = 0.238 ±0.005 (5) 

This calculation is shown in Fig. 4. In terms of a sin- 
gle parameter constraint on N v it is N v — 2. (95 
% conf.). Compared to the recent calculation by Abaza- 
jian [62| of a BBN-only constraint of 1.7 < N v < 3.5 
(95 % conf.) this is a significant improvement. Very 
interestingly the new limit suggests the possibility that 
N v is actually less than 3. This is for instance possible 
in scenarios with extremely low reheating temperature 

Of course this conclusion is mainly based on the fact 
that CMB and LSS data prefers a slightly higher value 
of than BBN. It should also be stressed that the 

estimates of the primordial abundances could be biased 
by systematic effects so that the quoted statistical error 
bar is not really meaningful. Therefore it is probably 
premature to talk of any inconsistency between the N v = 
3 prediction of the standard model and observations. 

In fact the argument can also be reversed. If N v is fixed 
to the standard model value of 3 then then CMB and LSS 
constraint on Vl^h 2 provides an accurate measurement of 
primordial He-4. Using the derived constraint on fli,h 2 
the 95% confidence range for Yp is 

0.2458 < Y P < 0.2471. (6) 

This could point to a serious underlying systematic effect 



in observational determinations of Yp, as discussed in 
Ref. 

Late entropy production - However, if entropy is pro- 
duced subsequent to BBN and prior to CMB formation 
it still behaves like an additional N v for CMB and LSS 
data, while having no effect on BBN. In the case where 
the decay of a massive particle produces the entropy the 
bound on N v can be translated into a bound on the mass 
and lifetime of the massive particle. The effective N v 
produced by complete decay of a non-relativistic species 
at lOkeV > T > 1 eV is 

N v , eS ~ 3 + 1.2 x (—) mt 7 v r y 2 / 3 , (7) 

where n Ut o is the number density of a standard model 
neutrino. 

The bound on N v therefore translates into a bound of 

* 3.3 (8) 
for any hypothetical massive particle. 

C. Neutrino lepton asymmetry and N v 

Whether the extra relativistic energy density is in the 
form of a non-zero chemical potential has no influence on 
CMB or LSS calculations since they are not sensitive to 
the flavour content of the additional energy density. How- 
ever, from a BBN perspective neutrinos in the electron 
flavour are very different from muon and tau neutrinos. 
The reason is that electron neutrinos and antineutrinos 
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enter directly into the weak interactions which maintain 
neutron-proton equilibrium. Adding electron neutrinos 
shifts the balance in the direction of a lower n/p ratio, 
i.e. the opposite effect of increasing N v . Therefore a large 
chemical potential in the muon or tau neutrino flavours 
can be hidden by compensating its effect on BBN with 
a small electron neutrino chemical potential. In fact, the 
additional N v would not necessarily have to be in the 
form of a non-zero or v T chemical potential. A large 
N v from additional weakly interacting species could also 
be compensated by a small electron neutrino chemical 
potential. 

However, since CMB and LSS observations are insen- 
sitive to flavour this degeneracy can be broken by com- 
bining BBN, CMB and SNI-a data ,41]. The effect of 
doing this can be seen in Fig. 5. From the present data a 
bound on the chemical potential in non-electron flavour 
neutrinos is |r/| < 2.3 at (95% conf.), where rj = fJt/T. 
In terms of N v the bound is N u < 6.5. As can be seen 
from the figure this bound comes almost entirely from 
the CMB+LSS data alone. 

It should be noted that any degree of neutrino mixing 
will act to equilibrate the flavours and possibly eliminate 
this masking effect. For small mass differences neutrino 
oscillations occur after weak freeze-out and consequently 
there is no effect on BBN. However, given the mass dif- 
ferences and close to maximal mixing angles suggested 
by the solar and atmospheric neutrino data, all three 
flavours will be almost equilibrated before weak freeze- 
out. This means that the bound on the lepton asymme- 
try in all flavours will be close to the BBN bound for 
the electron sector, i.e. < 0.15 EH EH [73, E3 
Therefore the bound derived here is not competitive, but 
is of course very robust because it does not depend on a 
knowledge of mass differences or mixing angles. It is also 
more robust in the sense that a large N v in hypothetical 
new particles can still be invisible to BBN if all neutrino 
species have a small chemical potential. Therefore the 
N u bound derived in this section does not have to be in 
the form of a f MjT chemical potential. 

D. Combining ^2 m » with N v - constraining LSND 

From the analyses in the above two sections it was 
found that: (a) An increasing ^ m v can be compensated 
by a decreasing Ho and (b) An increasing N v can be 
compensated by an increasing Hq. One might therefore 
wonder whether a model with non-zero m„, combined 
with N v > 3 can provide a good fit to the data. 

In order to test this we have performed a likelihood 
analysis for ^ m„ for different values of N v . We show 
this in Fig. 6. This analysis was performed with all avail- 
able data and priors. 

As can be seen from the figure, the best fit actually is 
actually shifted to higher ^ m„ when N v increases, and 
the conclusion is that a model with high neutrino mass 
and additional relativistic energy density can provide ac- 



10 




0.020 0.025 0.030 



FIG. 5: The full lines show 68% and 95% confidence regions 
in the (Qbh 2 , N v ) plane for the case where the additional N v 
is compensated during BBN by a small u e chemical potential. 
The full contours are for BBN data alone, whereas the dashed 
lines are for CMB and LSS data. 















/ 


/- 

/ - 

/ : 






/ 

/ 






/ . ■ ' ^ 















0.0050.010.0150.020.0250.03 

n^ 2 



FIG. 6: Ax 2 as a function of ^ m„ for various different values 
of N v . The full line is for iV„ = 3, the dotted for N v = 4, and 
the dashed for N v = 5. A\ 2 is calculated relative to the best 
fit N v — 3 model. 

ceptable fits to the data. As a function of N v the upper 
bound on ^2 m v is (at 95% confidence) 

X)m v <1.01eV for7V„ = 3 
X)m„<1.38eV for N v = 4 
£m„< 2.12 eV for N v = 5 

This has significant implications for attempts to con- 
strain the LSND experiment using the present cosmo- 
logical data. Pierce and Murayama conclude from the 
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present MAP limit that the LSND result is excluded [85| 
(see also Ref. |8(j). 

However, for several reasons this conclusion does not 
follow trivially from the present data. In general the 
three mass differences implied by Solar, atmospheric 
and the LSND neutrino measurements can be arranged 
into either 2+2 or 3+1 schemes. Recent analyses [93 
of experimental data have shown that the 2+2 mod- 
els are ruled out. The 3+1 scheme with a single mas- 
sive state, 7714, which makes up the LSND mass gap, is 
still marginally allowed in a few small windows in the 
(Am 2 , sin 2 29) plane. These gaps are at (Am 2 , sin 2 29) ~ 
(0.8 eV 2 , 2 x 10~ 3 ), (1.8 eV 2 , 8 x 10~ 4 ), (6eV 2 , 1.5 x 10~ 3 ) 
and (10eV 2 , 1.5xl0 -3 ). These four windows corresponds 
to masses of 0.9,1.4,2.5 and 3.2 eV respectively. From 
the Solar and atmospheric neutrino results the three light 
mass eigenstates contribute only about 0.1 eV of mass if 
they are hierarchical. This means that the sum of all 
mass eigenstate is close to 7714. 

The limit for N u = 4 which corresponds roughly to the 
LSND scenario is J2 m v ~ 1-4 eV, which still leaves the 
lowest of the remaining windows. The second window at 
to ~ 1.8 eV is disfavoured by the data, but not at very 
high significance. 

The second reason why the LSND sterile neutrino is 
not necessarily ruled out is that it was not necessarily 
fully equilibrated in the early universe. The cosmological 
mass limit applies only to a species with decoupling tem- 
perature around 1 MeV. If a sterile species is not fully 
equilibrated its number density is lower than that of a 
standard active neutrino species and therefore the mass 
limit will be less restrictive. In the simplest possible case 
a neutrino in a the lowest of the four LSND windows 
will have an abundance of roughly 0.5 times that of a 
standard active neutrino, whereas neutrino masses and 
mixings corresponding to the other three windows will 
lead to almost perfect equilibration [87l l8Sj. However, 
for instance the presence of a non-zero lepton number 
can strongly suppress the production of sterile neutrinos 

IV. DISCUSSION 

We have calculated improved constraints on neutrino 
masses and the cosmological relativistic energy density, 
using the new WMAP data together with data from the 
2dFGRS galaxy survey. 

Using CMB and LSS data together with a prior from 
the HST key project on Hg yielded an upper bound of 
J2 m i> < 1-01 eV (95% conf.). While this excludes most 



of the parameter range suggested by the claimed evidence 
for neutrinoless double beta decay in the Heidelberg- 
Moscow experiment, it seems premature to rule out this 
claim based on cosmological observations. 

Another issue where the cosmological upper bound on 
neutrino masses is very important is for the prospects of 
directly measuring neutrino masses in tritium endpoint 
measurements. The successor to the Mainz experiment, 
KATRIN, is designed to measure an electron neutrino 
mass of roughly 0.25 eV [9l|. or in terms of the sum of 
neutrino mass eigenstates, ^ to„ < 0.75 eV. The WMAP 
result of m v ~ 0-7 e V (95% conf.) already seems to ex- 
clude a positive measurement of mass in KATRIN. How- 
ever, this very tight limit depends on priors, as well as 
Ly-a forest data, and the more conservative present limit 
of ^2, m v ~ 1.01 eV (95% conf.) does not exclude that 
KATRIN will detect a neutrino mass. 

From the data we also inferred a limit on N v of N u = 
4-0^2 1 (95% conf.) on the equivalent number of neutrino 
species. This is a marked improvement over the previous 
best limit of roughly N„ < 13 [io|. 

When light element measurements of He-4 and D are 
included the bound is strengthened considerably to N v = 
2.6^o!3 (95 % conf.). Interestingly this suggests a possible 
value of N v which is less than 3. This could be the case 
for instance in scenarios with very low reheating tem- 
perature where neutrinos were never fully equilibrated 
HE EE E3- However, it should be stressed that pri- 
mordial abundances could be dominated by systematics. 
Therefore it is probably premature to talk of a new BBN 
"crisis" . 

The bound on N v also translates into a bound on a 
possible neutrino lepton asymmetry in non-electron neu- 
trinos of \r) v \ < 2.4. Even though this bound is much 
stronger than the previous bound of 2.6 [Il| it is weak 
compared to the bound obtainable from BBN consider- 
ations alone when flavour oscillations are accounted for 

El IE IE EE E3 

Finally, we also found that the neutrino mass bound 
depends on the total number of light neutrino species. 
In scenarios with sterile neutrinos this is an important 
factor. For instance in 3+1 models the mass bound in- 
creases from 1.0 eV to 1.4 eV, meaning that the LSND 
result is not ruled out by cosmological observations yet. 
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